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ABSTRACT 


<^?The  structural  dynamic  response  of  turbomachinery  components 
to  aero-thermodynamic  distortion  induced  excitations  is  iw—tberiF— - 
of  major  concern  in  the  design  of  advanced  gas  turbine  engines. 
The  rotor  speeds  at  which  these  resonant  forced  responses  occur 
can  be  predicted  with  Campbell  diagrams.  However,  due  to  the 


inadequacies  of  the  existing  time-variant  aerodynamic  models,  no 
accurate  prediction  can  currently  be  made  for  the  amplitude  of 
the  resulting  stresses. 


The  overall  objective  of  this  research  program  is  to  quanti¬ 
tatively  investigate  the  fundamental  phenomena  relevant  to  aero- 
thermodynamic  distortion  induced  structural  dynamic  blade 
response  in  multi-stage  gas  turbine  fans  and  compressors  Unique 
unsteady  aerodynamic  data  will  be  obtained  to  validate  and  indi¬ 
cate  necessary  refinements  to  state-of-the-art  analyses  and  to 
direct  the  modeling  of  new  analyses.  Also,  for  the  first  time,>a 
first  principles  capability  to  predict  the  vibrational  response 
amplitude  of  blading  due  to  aerodynamic  excitations  will  be 
developed.  Accession  FoT  | 


Ia-i 


This  report  summarizes  the  progress  and  results  obtained 
during  the  period  April  16,  1983  to  April  15,  1984.  These 
include.:,  the  calibration  of  stator  vanes  instrumented  with 
dynamic  pressure  transducers;  the  completion  and  verification  of 
the  steady  and  dynamic  data  acquisition  and  analysis  procedures; 
the  check-out  and  initial  experimental  study  of  the  first  stage 
vane  row  unsteady  aerodynamics ;  the  initiation  of  the  development 
of  the  dynamic  instrumentation  and  calibration  procedures  for  the 
rotor  studies;  and  the  completion  of  the  structural  dynamics  por¬ 
tion  of  the  coupled  mode  blade  vibration  energy  balance  struc¬ 
tural  dynamics  math  model,  f  -  '  -  7  '  ' ' "  '''  / 
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INTRODUCTION 


I . 


The  structural  dynamic  response  of  fern,  compressor,  and  tur¬ 
bine  blading  to  aero-thermodynamic  distortion  induced  excitations 
is  an  item  of  rapidly  increasing  concern  to  designers  and 
manufacturers  of  gas  turbine  engines  for  advanced  technology 
applications.  Destructive  aerodynamic  forced  responses  of  gas 
turbine  engine  blading  have  been  generated  by  a  wide  variety  of 
aero-thermodynamic  distortion  sources.  These  include:  blade 
wakes;  multi-staging  interaction  effects;  large  angle  of  attack 
or  yaw;  engine  exhaust  recirculation;  cross-flow  at  the  inlet; 
pressure  variations  on  the  engine  due  to  external  aerodynamics; 
and  armament  firing.  The  following  physical  phenomena  have  been 
identified  as  being  significant  with  regard  to  aero-thermodynamic 
distortion  induced  structural  dynamic  responses  of  fan  and 
compressor  blading. 

*  Resonance 

*  Multi-stage  interactions 

*  Stall 

*  Inlet  gusts 


Time-varying  inlet  flows 


*  Plow  separation 


*  Potential  effects 

*  Turbulence 


The  first  principles  prediction  of  the  structural  dynamic 
response  associated  with  all  of  the  above  physical  phenomena  is 
identical  and  involves  the  following  elements.  Spatially 
periodic  variations  in  pressure,  velocity,  temperature,  and  flow 
direction  of  the  exit  flow  field  of  an  upstream  element  appear  as 
temporally  varying  in  a  coordinate  system  fixed  to  the  downstream 
row.  As  a  result,  individual  airfoils  are  subject  to  a  time- 
variant  aerodynamic  forcing  function  which  can  induce  high  level 
vibratory  stresses. 

The  analysis  of  the  aerodynamically  forced  response  vibra¬ 
tory  behavior  of  a  blade  or  vane  row  requires  a  definition  of  the 
unsteady  forcing  function  in  terms  of  its  harmonics.  The  time- 
variant  aerodynamic  response  of  the  airfoil  to  each  harmonic  of 
this  forcing  function  is  then  assumed  to  be  comprised  of  two 
parts.  One  is  due  to  the  disturbance  being  swept  past  the  non¬ 
responding  fixed  airfoils.  The  second  arises  when  the  airfoils 
respond  to  this  disturbance.  Mathematically  these  effects  are 
modeled  by  two  distinct  analyses.  A  linearized  small  perturba¬ 
tion  gust  analysis  is  used  to  predict  the  time-variant  aerodynam¬ 
ics  of  the  fixed  non-responding  airfoils  to  each  harmonic  of  the 
disturbance.  A  self-induced  unsteady  aerodynamic  analysis 


wherein  the  airfoils  are  assumed  to  be  harmonically  oscillating 
is  then  used  to  predict  the  additional  aerodynamic  effects  due  to 
the  airfoil  response.  Superposition  of  these  two  effects  can  be 
performed  only  with  knowledge  of  the  modal  pattern  and  amplitude 
of  response  of  the  blading  because  the  magnitude  of  the  pressure 
field  resulting  from  the  airfoil  motion  is  dependent  upon  the 
amplitude  of  this  motion.  Thus  a  solution  with  key  elements  con¬ 
sisting  of  the  gust  analysis,  a  self -induced  unsteady  aerodynamic 
analysis,  and  an  airfoil  structural  analysis  is  necessary  to 
predict  the  total  response  of  an  airfoil  subjected  to  an  upstream 
generated  spatially  periodic  disturbance. 

The  unsteady  small  perturbation  gust  and  self- induced 
unsteady  aerodynamic  analyses  are  two-dimensional  and,  as  such, 
are  coupled  to  the  airfoil  structural  analysis  by  means  of  a 
strip  theory  approximation.  Thus,  the  airfoil  is  considered  to 
consist  of  a  series  of  individual  and  independent  two-dimensional 
aerodynamic  regions.  The  time-variant  aerodynamic  analyses  are 
then  applied  to  each  such  individual  region,  with  the  charac¬ 
teristic  parameters  including  the  Mach  number,  reduced  frequency, 
stagger  angle,  and  solidity,  taken  as  the  average  value  at  the 
inlet  boundary.  It  should  be  noted  that  there  is  no  coupling 
between  adjacent  aerodynamic  regions.  Hence,  aerodynamic  forced 
vibrations  involving  spanwise  variations  in  unsteady  aerodynamics 
cannot  be  treated,  i.e.,  if  spanwise  variations  exist  in  the 
aerodynamic  forcing  function  not  caused  by  simple  inlet  spanwise 
variations  in  the  Mach  number,  reduced  frequency,  stagger  angle, 


solidity,  such  a  strip  theory  design  system  is  of  no  value, 
of  the  previously  noted  forced  response  sources  may  fit  into 
category,  as  for  example,  the  case  of  rotor  tip  vortices 
ating  a  forced  response  in  a  downstream  blade  or  vane  row. 


II.  PROGRAM  OBJECTIVES 


Overall  Objectives 


The  overall  objective  of  this  research  program  is  to  quanti¬ 
tatively  investigate  the  fundamental  phenomena  relevant  to  dis¬ 
tortion  generated  aero-thermodynamic  induced  structural  dynamic 
effects  in  gas  turbine  engine  blade  rows.  Unique  unsteady  aero¬ 
dynamic  data  will  be  obtained  on  stationary  and  rotating  blade 
rows  to  validate  and  indicate  necessary  refinements  to  current 
state-of-the-art  analyses  and  to  direct  the  modeling  and  develop¬ 
ment  of  new  analyses.  Also,  for  the  first  time,  a  first  princi¬ 
ples  capability  to  predict  the  vibrational  response  amplitude  of 
blading  due  to  aerodynamic  excitations  based  on  an  energy  balance 
technique  will  be  developed. 

From  first  principles  considerations,  the  relevant  fundamen¬ 
tal  physical  phenomena  are  identical  for  the  various  sources  of 
aero-thermodynamic  distortion.  Hence,  to  accomplish  this  overall 
objective  in  a  timely  and  efficient  ma-mer ,  while  obtaining 
results  of  direct  interest  and  significance  to  the  gas  turbine 
engine  community,  this  proposed  research  program  is  concerned 
with  the  time-variant  aerodynamics  and  structural  dynamic 
response  of  multi-stage  stationary  and  rotating  blade  rows,  with 
the  primary  source  of  excitation  initially  being  the  wakes  from 
upstream  blade  elements.  The  overall  objective  includes  the  fol- 
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The  experimental  determination  of  the 
fundamental  time-variant  gust  aerodynamics 
associated  with  variations  in  incidence 
angle  (loading),  reduced  frequency, 
solidity,  and  multi-stage  effects  on 
both  stationary  and  rotating  blade  rows 
as  well  as  the  investigation  of  the 
validity  of: 

*  the  two-dimensional  linearity  and 
superpositioning  assumptions; 

*  the  small  perturbation  modeling  concept; 

The  development  of  a  first  principles 
state-of-the-art  capability  to  predict  the 
vibrational  response  amplitude  of  blading 
based  on  an  energy  balance  technique,  with 
key  elements  consisting  of  the  following 
coupled  existing  analyses: 

*  a  flat  plate  cascade  transverse  gust 
analysis ; 

*  a  flat  plate  cascade  self-induced 
unsteady  aerodynamics  analysis; 

*  an  airfoil  structural  analysis; 


*  The  development  of  a  subsonic  compressible 
unsteady  aerodynamic  cascade  analysis. 

Thus,  this  program  is  directed  at  providing  fundamental 
time-variant  aerodynamic  data  whic''  will  not  only  validate  and 
indicate  refinements  to  the  current  state-of-the-art  two- 
dimensional  gust  analyses  but  will  also  address  the  validity  of 
the  most  basic  assumptions  inherent  in  these  analyses  and  in  the 
structure  of  forced  response  design  systems.  In  addition  a 
state-of-the-art  first  prinicples  predictive  aerodynamically 
forced  response  analysis  based  on  an  energy  balance  technique 
will  be  developed. 

Specific  Objectives 

The  following  were  the  specific  objectives  for  thi3  report¬ 
ing  period. 


Stationary  Vane  Row  Experiments 

*  Calibration  of  the  stator  vanes  instrumented 
with  Kulite  dynamic  pressure  transducers 

*  Completion  and  verification  of  the  steady  and 
dynamic  data  acquisition  and  analysis  procedures 

*  On-line  check  out  of  the  data  acquisition 
system 


*  Initiation  of  the  acquisition  and  analysis  of 
the  first  stage  vane  row  time-variant  aerodynami 
data 

Rotor  Blade  Studies 

*  Initiation  of  the  development  of  the  dynamic 
instrumentation  and  calibration  procedures  for 
the  first  stage  rotor  blade  row 

Analysis 


I 


i 


Parametric  studies  completed  with  the  coupled 
mode  blade  vibration  energy  balance  structural 
dynamics  mathematical  model 


III.  TECHNICAL  APPROACH 


The  technical  approach  to  achieve  the  overall  program  objec¬ 
tives  requires  that  extensive,  high-quality,  detailed  aerodynamic 
data  be  acquired  from  benchmark  experiments  which  model  the  fun¬ 
damental  flow  physics  of  aero-thermodynamic  induced  structural 
dynamic  effects  in  gas  turbine  engines.  These  data  must  be 
analyzed  and  correlated  with  appropriate  state-of-the-art  ana¬ 
lyses  to  ascertain  their  range  of  validity  and/or  to  indicate 
necessary  refinements.  In  addition,  advanced  mathematical  models 
and  techniques  for  the  prediction  of  these  phenomena  will  be 
developed  and  experimentally  verified. 

The  approach  to  achieving  these  experimental  objectives  is 
to  measure  the  time-variant  pressure  distribution  in  controlled 
benchmark  experiments  which  model  the  fundamental  flow  physics 
and  thereby  identify  and  quantify  the  key  unsteady  aerodynamic 
parameters  relevant  to  aero-thermodynamic  distortion  induced 
response  of  blading.  In  particular,  this  research  program  is 
aimed  at  measuring  the  fundamental  unsteady  gust  aerodynamics 
associated  with  variations  in  incidence  angle  (loading),  reduced 
frequency,  solidity,  multi-stage  effects,  as  well  as  investigat¬ 
ing  the  validity  of  the  two-dimensional  linearity  and  superposi¬ 
tioning  assumptions  and  the  applicability  of  the  small  perturba¬ 
tion  modeling  concept.  This  requires  obtaining  airfoil  surface 
chordwise  and  spanwise  unsteady  pressure  data  on  both  stationary 


and  rotating  blade  rows  as  well  as  detailed  inlet  and  exit  plane 
flow  surveys  in  a  multi-stage  axial  flow  compressor.  These  data 


will  be  correlated  with  state-of-the-art  unsteady  cascade 
transverse  gust  analyses  and  will  be  used  to  determine  their 
range  of  validity  and  to  direct  the  development  of  an  advanced 
first  principles  forced  vibration  aerodynamic  analysis  which  the 
designer  can  use  with  confidence.  In  addition,  a  unique  first 
principles  predictive  aerodynamically  forced  response  analysis 
will  be  developed  based  on  an  energy  balance  technique. 


Facility 

In  the  first  stage  of  a  multi-stage  axial  flow  compressor 
the  wakes  from  the  upstream  rotor  blades  are  the  source  of  the 
aerodynamically  induced  fluctuating  surface  pressures  on  the 
downstream  stator  vanes;  i.e.,  the  rotor  blade  wakes  define  the 
forcing  function  to  the  downstream  stator  vane  row.  In  the  aft 
stages  of  the  compressor  the  wakes  from  all  of  the  upstream  blad¬ 
ing,  both  rotors  and  stators,  are  the  source  of  the  aerodynami¬ 
cally  induced  fluctuating  surface  pressure  distributions  on  the 
downstream  stator  vane  and  rotor  blade  rows.  Hence,  it  is  neces¬ 
sary  to  experimentally  model  all  of  the  significant  features 
which  define  this  forcing  function.  These  include  the  variation 
of  incidence,  the  wave  form,  the  velocity  (pressure)  variation, 
the  various  blade  row  interactions  in  a  multi-stage  machine,  and 
the  reduced  frequency  (k  -  wC/  2V  .  .)  which  defines  the  order 

aX  1  cL  X 

of  magnitude  of  the  unsteadiness.  (Forced  response  problems  in 
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turbomachines  typically  have  first  harmonic  reduced  frequency 
values  varying  from  0.2  to  8.0.)  The  above  noted  features  can 
all  be  simulated  in  the  Purdue  University  Multi-Stage  Axial  Flow 
Research  Compressor. 

Figure  1  presents  a  schematic  of  the  overall  facility,  with 
the  compressor  test  section  shown  in  greater  detail  in  Figures  2 
and  3.  The  compressor  is  driven  by  a  15  HP  DC  electric  motor 
with  a  speed  range  of  300  to  3000  RPM.  The  inlet  section  is 
located  aft  of  the  drive  motor.  In  the  exit  of  this  section  are 
located  38  variable  geometry  guide  vanes  which  direct  the  flow 
into  the  test  section.  Three  identical  compressor  3tages  are 
mounted  in  the  test  section,  which  has  an  annulus  with  constant 
hub  (0.300  m)  and  tip  (0.420  m)  diameters.  The  exit  flow  from 
the  test  section  is  directed  through  a  series  of  flow  straight- 
eners  into  a  venturi  meter  which  enables  the  mass  flow  rate  to  be 
determined.  To  throttle  the  compressor,  an  adjustable  plate  is 
located  at  the  exit  of  the  diffuser  of  the  venturi. 

Each  of  the  three  identical  compressor  stages  consists  of  43 
rotor  blades  and  41  stator  vanes.  Hence,  the  interblade  phase 
angle  for  these  experiments  is  17.56°.  These  free  vortex  design 
airfoils  have  a  British  C4  section  profile,  a  chord  of  30  mm,  and 
a  maximum  thickness-to-chord  ratio  of  0.10.  The  overall  airfoil 
and  compressor  characteristics  are  presented  in  Table  I. 


gure  3.  Close-up  view  of  the  blading  of  the 
Three-Stage  Research  Compressor 
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Type  of  Airfoil 
Number  of  Blades 
Chord,  C(ram) 

Solidity,  C/S 
Camber 

Aspect  Ratio 
Thickness/Chord  (%) 

Flow  Rate  (kg/second) 

Rotor-Stator  Axial  Spacing  (mm) 
Design  Axial  Velocity  (meters/sec) 
Rotational  Speed  (RPM) 

Number  of  Stages 
Stage  Pressure  Ratio 
Inlet  Tip  Diameter  (mm) 

Hub/Tip  Radius  Ratio 
Stage  Efficiency  (%) 

TABLE  I.  AIRFOIL  MEAN  SECTION 


Rotor 

Stator 

C4 

C4 

43 

41 

30 

30 

1.14 

1.09 

27.95 

27.70 

2.0 

2.0 

10 

10 

2.66 

22.2 

32.0 

3000 

3 

1.003 

420 

0.714 

85 

CHARACTERISTICS  AND 


COMPRESSOR  DESIGN  POINT  CONDITIONS 


Conventional  steady-state  instrumentation  is  used  to  deter¬ 
mine  the  flow  properties  in  the  compressor.  The  inlet  tempera¬ 
ture  is  measured  by  four  equally  spaced  thermocouples  at  the 
inlet  of  the  compressor.  Casing  static  taps,  equally  spaced  cir¬ 
cumferentially,  allow  the  measurement  of  the  static  pressure 
between  each  blade  row.  Traversing  gear  instrument  stations  pro¬ 
vided  between  each  blade  row  are  used  to  measure  the  mean  flow 
incidence  angle  on  the  first-stage  stator  vanes .  A  thermocouple 
and  Kiel  probe  located  downstream  are  used  to  measure  the 
compressor  exit  temperature  and  total  pressure,  respectively. 
The  mass  flow  is  measured  with  the  calibrated  venturi  meter 
located  downstream  of  the  compressor  test  section.  A  shaft 
mounted  60  tooth  gear  and  a  magnetic  pickup  provide  the  rotor 


IV.  PROGRAM  STATUS  AND  RESULTS 


All  of  the  previously  noted  specific  objectives  established 
for  this  reporting  period  have  been  successfully  met. 


Stationary  Vane  Row  Experiments 

The  calibration,  steady  and  dynamic  data  acquisition  and 
analysis  procedures,  and  the  check-out  of  the  data  acquisition 
system  have  been  completed.  Also,  the  initial  acquisition  and 
analysis  of  the  first  stage  vane  row  time-variant  data  has  begun. 


Initial  Check-Out  Experiment 

In  this  initial  first  stage  stator  vane  check-out  study,  the 
wakes  from  the  first  stage  rotor  blades  are  the  source  of  the 
aer odynamically  induced  fluctuating  surface  pressure  distribu¬ 
tions,  i.e.,  the  rotor  wakes  define  the  forcing  function  to  the 
first  stage  stator  vanes,  as  schematically  depicted  in  Figure  4. 

The  unsteady  data  of  fundamental  interest  are  the  chordwise 
distribution  of  the  complex  time-variant  pressure  difference 
across  the  chordline  of  the  first  stage  stator  vane.  This 
requires  the  measurement  of  the  dynamic  pressure  distributions  on 
the  pressure  and  suction  surfaces  of  the  first  stage  stator  vane 


Figure  4.  Schematic  of  flow  field 
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at  identical  chordwise  locations-  These  airfoil  surface  fluc¬ 
tuating  pressure  measurements  are  accomplished  by  instrumenting 
pairs  of  stator  vanes  with  Kulite  thin-line  design  dynamic  pres¬ 
sure  transducers. 

To  minimize  any  flow  disturbances  generated  by  the  transduc¬ 
ers,  they  were  embedded  in  the  airfoil  and  connected  to  the  sur¬ 
face  by  a  static  pressure  tap.  The  transducer  lead  wires  were 
placed  in  milled  slots  and  carried  out  through  hollow  trunnions. 
A  similar  installation  was  used  for  the  transducers  on  the  pres¬ 
sure  surface  at  the  same  chord  locations.  Thus,  in  this  initial 

H. 

check-out  study,  sets  of  two  instrumented  vanes  were  used,  one 
with  suction  surface  transducers  and  one  with  pressure  surface 
transducers.  These  vanes  were  located  such  that  a  complete  flow 
passage  was  instrumented.  Conventional  static  pressure  taps  were 
also  provided  at  the  same  locations.  These  static  pressure  taps 
provide  the  steady  state  aerodynamic  load  distribution  on  the 
airfoil. 

Both  steady-state  and  time-variant  data  are  being  acquired. 
Figure  5  presents  a  schematic  of  the  complete  data  acquisition 
and  analysis  system  which  is  centered  on  a  PDP  11-23  computer. 

The  steady-state  data  of  interest  include:  (1)  the  static 
pressures  on  the  stator  vane  surfaces,  accomplished  with  chord- 
wise  distributions  of  surface  static  pressure  taps;  (2)  the  pres¬ 
sures  and  temperatures  throughout  the  compressor,  thereby  permit¬ 
ting  the  compressor  operating  map  to  be  determined,  as  shown  in 
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Pigure  6.  The  compressor  map  is  used  to  define  the  operating 
points,  in  terms  of  overall  pressure  ratio  and  corrected  mass 
flow  rate,  at  which  the  unsteady  surface  pressure  measurements 
are  obtained. 

Steady-state  data  acquisition  follows  the  standard  evalua¬ 
tion  procedure.  At  the  selected  corrected  speed,  the  compressor 
is  stabilized  for  approximately  10  minutes,  after  which  the 
steady-state  sequence  is  initiated.  The  data  are  then  scanned, 
and  the  corrected  mass  flow,  pressure  ratio,  and  corrected  speed 
determined . 

The  time-variant  data  acquisition  and  analysis  technique 
used  is  based  on  a  data  averaging  or  signal  enhancement  concept. 
The  key  to  such  a  technique  is  the  ability  to  sample  data  at  a 
preset  time.  In  this  investigation,  the  data  of  interest  are 
being  generated  at  the  blade  passing  frequency.  Hence,  an  opti¬ 
cal  encoder,  delivering  a  square  wave  voltage  signal  with  a  dura¬ 
tion  in  the  microsecond  range,  was  mounted  on  the  rotor  shaft  and 
used  as  the  time  or  data  initiation  reference  to  trigger  the  A-D 
multiplexer  system.  This  system  is  capable  of  digitizing  signals 
simultaneously  at  rates  to  5  megahertz  per  channel,  storing  2048 
data  points  per  channel. 

In  this  initial  check-out  experiment,  the  effect  of  averag¬ 
ing  the  time-variant  pressure  signals  from  the  blade  mounted 
dynamic  pressure  transducers  was  considered.  Figures  7  through 
11  display  the  time-variant  pressure  signal  from  the  14.1%  chord 
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Figure  6.  Overall  compressor  performance  map 
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pressure  signal  averaged  over  1  rotor  revolution 
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averaged  pressure  signal  averaged  over  25  rotor  revolutions 
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pressure  signal  averaged  over  50  rotor  revolutions 
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pressure  surface  dynamic  pressure  transducer  for  1  rotor  revolu¬ 
tion  and  averaged  over  25,  50,  75,  and  100  rotor  revolutions, 
respectively.  As  seen,  averaging  greatly  reduces  the  random 
fluctuations  superimposed  on  the  harmonic  pressure  signal  of 
interest.  Also,  these  time-variant  pressure  signals  are  essen¬ 
tially  unchanged  when  averaged  over  75  or  more  rotor  revolutions, 
Figures  10  and  11.  Thus,  signal  averaging  for  100  rotor  revolu¬ 
tions  are  used  for  the  time-variant  data  acquisition  of  this  ini¬ 
tial  check-out  experiment. 

At  each  steady-state  operating  point,  an  averaged  time- 
variant  data  set,  consisting  of  the  Kulite  dynamic  pressure 
transducer  signals  digitized  at  a  rate  of  20  KHz  and  averaged 
over  100  rotor  revolutions,  was  obtained.  These  rotor  revolu¬ 
tions  were  not  consecutive  due  to  the  finite  time  required  for 
the  A-D  multiplexer  system  to  sample  the  data  and  the  computer  to 
then  read  the  digitized  data.  Figures  12  and  13  present  typical 
examples  of  the  digitized  dynamic  pressure  signals  at  14.1%  chord 
on  the  vane  pressure  and  suction  surfaces,  respectively. 

Each  of  these  digitized  signals  was  Fourier  decomposed  into 
harmonics  by  means  of  an  FFT  algorithm.  Figure  14  shows  an  exam¬ 
ple  of  this  decompos it ion  for  the  14.1%  chord  pressure  surface 
transducer  signal  at  73%  span.  As  seen,  the  transducer  signal 
contains  a  dominant  fundamental  frequency  equal  to  the  blade 
passing  frequency  and  a  smaller  second  harmonic  component.  In 
addition,  the  averaged  signal  exhibits  minimal  non-harmonic  con¬ 
tent.  From  this  Fourier  decomposition,  both  the  magnitude  of 
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Fourier  analysis  of  averaged  pressure  signal  at 
14.1%  chord  on  the  pressure  surface 


each  component  and  its  phase  lag  as  referenced  to  the  optical 
encoder  pulse  can  be  determined. 


Cascade  gust  unsteady  aerodynamic  analyses  predict  the 
unsteady  pressure  difference  across  the  vane  as  a  function  of 
airfoil  chord  in  the  form  of  a  dynamic  pressure  difference  coef¬ 
ficient  and  an  aerodynamic  phase  lag.  The  dynamic  pressure  coef¬ 
ficient  is  normalized  with  respect  to  the  steady-state  properties 
of  the  flow  and  the  magnitude  of  the  transverse  gust.  The  aero¬ 
dynamic  phase  lag  is  referenced  to  a  transverse  gust  at  the  lead¬ 
ing  edge  of  the  reference  airfoil.  Hence,  to  enhance  the 
experiment-theory  correlation  process,  the  data  were  adjusted  in 
phase  and  magnitude  to  agree  with  the  predictions  at  the  leading 
edge  transducer  location.  Thus,  the  final  form  of  the  dynamic 
data  consists  of  a  dynamic  pressure  coefficient  and  an  aero¬ 
dynamic  phase  lag.  It  should  be  noted  that  this  adjustment  in 
phase  and  magnitude  will  not  be  necessary  in  the  ensuing  experi¬ 
ments  as  a  crossed  hot-wire  system  will  be  used  to  quantify  the 
forcing  function. 

Surface  static  pressure  measurements  on  the  49.1%  and  the 
73.0%  span  streamlines  of  the  first  stage  stator  vane  quantify 
the  steady-state  aerodynamic  chordwise  loading  distributions. 
Figure  15  presents  these  surface  chordwise  steady-state  static 
pressure  distributions. 

As  seen,  at  small  negative  incidence  angle  values,  the  aero¬ 
dynamic  loading  distributions  differ  over  the  front  portion  of 


Incidence 


the  chord,  but  coincide  over  the  aft  portion  for  these  two  span- 
wise  streamlines.  As  the  incidence  angle  is  increased  to  small 
positive  values,  the  pressure  surface  static  pressure  distribu¬ 
tions  coincide,  but  the  suction  surface  distributions  do  not.  At 
large  positive  incidence  angle  values,  the  two  span  streamline 
static  pressure  distributions  are  very  nearly  equal  except  for  a 
slight  deviation  in  the  leading  edge  portion  of  the  suction  sur¬ 
face  . 

Pigures  16  through  21  present  the  stator  vane  surface 
dynamic  pressure  coefficient  and  aerodynamic  phase  lag  data  for 
the  first  harmonic  of  the  unsteady  pressure  difference  across  the 
vane  as  a  function  of  percent  airfoil  chord.  Also  included  in 
these  figures  are  the  predictions  obtained  from  a  compressible 
flow,  small  perturbation,  flat  plate  airfoil  cascade,  transverse 
gust  unsteady  aerodynamic  analysis.  In  this  unique  analysis,  the 
effect  of  a  specified  steady-state  separation  on  the  airfoil  suc¬ 
tion  surface  is  modeled  by  a  thin  slip  region  which  extends  from 
the  specified  separation  point  chordwise  location  to  aft  of  the 
blade  trailing  edge. 

Figures  16  and  17  present  the  small  negative  and  positive 
incidence  angle  dynamic  pressure  data  on  the  73.0%  span  stream¬ 
line.  As  seen,  good  correlation  exists  between  the  data  and  the 
unseparated  predictions  for  the  dynamic  pressure  coefficient. 
However,  the  aerodynamic  phase  lag  data  exhibit  a  somewhat  dif¬ 
ferent  chordwise  distribution  than  that  of  the  prediction.  In 
particular,  the  aerodynamic  phase  lag  first  increases  and  then 
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decreases  as  a  function  of  chord,  with  the  decrease  being  more 
pronounced  as  the  incidence  angle  is  increased.  This  chordwise 
variation  in  the  aerodynamic  phase  lag  data  is  attributed  to  the 
camber  of  the  stator  vane.  The  difference  between  the  small 
negative  and  positive  incidence  angle  phase  lag  distributions  is 
due  to  the  differences  in  the  steady-state  aerodynamic  loading, 
as  detailed  in  Figure  15. 

The  small  negative  and  positive  incidence  angle  results  on 
the  49.1%  span  streamline  are  presented  in  Figures  18  and  19, 
respectively.  As  seen,  at  the  small  negative  incidence  angle, 
the  aerodynamic  phase  lag  data-theory  correlation  is  very  similar 
to  the  corresponding  data  on  the  73.0%  streamline.  At  the  posi¬ 
tive  incidence  angle,  a  larger  decrease  in  the  phase  lag  value  at 
63.7%  chord  is  found.  The  dynamic  pressure  coefficient  correla¬ 
tion  at  49.1%  span  is  generally  good  at  both  incidence  angle 
values.  However,  the  29.1%  chord  data  point  is  noticeably  lower 
in  value  than  the  corresponding  prediction  at  the  negative 
incidence  angle  and  increased  in  value  over  the  prediction  at  the 
positive  incidence  angle  value.  These  variations  in  the  aero¬ 
dynamic  phase  lag  and  the  dynamic  pressure  coefficient  data  are 
due  to  the  differences  in  the  steady-state  aerodynamic  loading, 
as  presented  in  Figure  15. 

Figures  20  and  21  present  the  large  positive  incidence  angle 
dynamic  pressure  data-theory  correlation  on  the  49.1%  and  73.0% 
span  streamlines,  respectively.  Also  included  in  these  figures 
are  the  correponding  unseparated  and  separated  flow  predictions. 


Both  sets  of  data  exhibit  analogous  trends.  As  per  the  lower 
incidence  angle  results,  the  aer odynamic  phase  lag  data  first 


increase  and  then  decrease  relative  to  the  predictions.  Also, 
the  aft  region  separated  flow  prediction  indicates  little  effect 
on  the  predicted  phase  lag  at  this  large  reduced  frequency  value. 
The  dynamic  pressure  coefficient  chordwise  data  distribution 
exhibits  the  same  trends  as  the  separated  flow  predictions.  In 
particular,  the  dynamic  pressure  coefficient  data  and  predictions 
decrease  as  the  chord  position  is  increased  over  the  front  of  the 
airfoil.  However,  while  the  nonseparated  prediction  continues  to 
decrease  over  the  aft  part  of  the  airfoil  chord,  the  data  and  the 
separated  flow  predictions  show  increased  values,  with  the  higher 
incidence  angle  data  exhibiting  the  largest  increase,  as  seen  in 
Figure  21. 


Rotor  Studies 

The  development  of  the  dynamic  instrumentation  and  calibra¬ 
tion  procedures  for  the  first  stage  rotor  blade  studies  has  been 
initiated.  Modifications  to  the  rotor  shaft  to  route  lead  wires 
to  the  slip  ring  assembly  has  been  initiated.  Also,  a  slip  ring 
assembly  has  been  selected  and  ordered.  In  addition,  design  has 
been  initiated  on  the  rotating  signal  conditioning  system  as  well 
as  the  dynamic  transducer  installation  on  the  rotor  blades. 


the  parametric  study  has  been  completed  with  the  coupled 
mode  blade  vibration  energy  balance  structural  dynamics  math 
model.  These  results  are  presented  in  the  Appendix. 


APPENDIX 


THE  COUPLED  RESPONSE  OF 
TURBOMACHINERY  BLADING  TO 
AERODYNAMIC  EXCITATIONS 


AIAA  Journal  of  Aircraft 
Volume  21,  No.  4,  April  1984 


The  Coupled  Response  of  Turbomachinery  Blading 
to  Aerodynamic  Excitations 


Daniel  Hoyniak*  and  Sanford  Fleeterf 
Purdue  University,  West  Lafayette,  Indiana 


Km  cwergy  balance  lecbniqae  is  developed  which  predicts  the  coupled  bending-tonion  mode  aerodynamic 
forced  response  at  an  airfoil.  The  effects  of  l he  various  aerody  namic  parameters  are  then  considered  utilizing  a 
subsonic  compressible  flow/ flat  plate  cascade  gust  analysis.  The  increased  coupling  between  the  torsion  and 
translation  modes  as  the  natural  frequencies  approach  one  another  is  shown.  It  is  also  demonstrated  that  the 
coopted- response  amplitudes  increase  with:  11  decreased  structural  damping,  2)  increased  solidity,  stagger  angle, 
and  Mach  nambers.  3)  ialerblade  phase  angles  corresponding  to  forward  (raveling  waves,  and  4)  shilling  of  the 
elastic  asrs  location  aft. 


Nomenclature 

a  =  distance  of  elastic  axis  aft  of  mid-chord 

A  -  unsteady  lift  influence  coefficient 

b  =  airfoil  semichord  (C/2) 

B  =  unsteady  moment  influence  coefficient 

C  =  airfoil  chord 

C/S  =  cascade  solidity  (airfoil  chord/airfoil  spacing) 
C,  =trf„  +  (A,.)* 

C,  =pjr.  +  (B»)" 

C,  =(p -m(oi,/u>-'  +  M,)"! 

C,  =(/l),),-p(u>,/i)-«» 

c,  =pr;n  -  (<*>„/*)']  +  (8.,)" 
c»  =  (B,,)' -vai(u>„/u>)2g„ 

C ,  =~{WC,U){AC)K 

c ,  =-(w0/unA0)' 

C,  =  -(WclU)(Bc)K 

c,„  =  ~{W0/U){AC)' 

D,  =  I  (C,C,  +C,Cj)  -  (C..M„),  +  C,  M„)'  1 1 
=|  (C.,C,-C4C»)-  |C.C)-(d„)/(fl,),|  I 

h  =  complex  translational  displacement 

h.,  -  amplitude  of  translational  oscillation 

i  =  v  -  I 

/ ,  =  mass  moment  of  inertia  per  unit  span  about  the 

elastic  axis 

it  =  reduced  frequency  (it  =  ub/  U) 

L  =  unsteady  lift 

in  =  mass  per  unit  span  of  the  airfoil 

1/  =  unsteady  moment 

•V„  =  |  (C,C,  +C»C,)  -  ICvM  , )'  +  C,C,„|  I 

V„„  =1  (C,C3-C,C,)-  [C,C\-C,„(A„)'|  | 

V =  I  ( C,„Cj  +  C.C, )  -  ( C  <  Bh ) 1  +  CjC, }  | 

v„,  =  i  ( c,c,  -  c4cl0 )  - 1 c-c,  - c,  ( a* ) '  1 1 

Re  [  |  =  real  part  of  (  I 

r ,  =  radius  of  gyration  about  elastic  axis 

S„  =  airfoil  static  moment  per  unit  span  about  the 
elastic  axis 
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U  =  freestream  velocity 

W  =  complex  transverse  gust  function 

Wc  =  amplitude  of  transverse  gust 

WORK  =  unsteady  work  per  cycle  of  oscillation 
x„  =  location  of  center  of  gravity  relative  to  elastic  axis 
a  =  complex  torsional  displacement 

a„  -  amplitude  of  torsional  motion 

n  =  mass  parameter 

p  -  fluid  density 

a  =  interblade  phase  angle 

=  airfoil  natural  frequency 
i  =  frequency  of  the  transverse  gust 

Subscripts 

G  =  gusi 

h  =  translation 

SI  =  self-induced  (aerodynamic  damping) 

u  =  torsion 

Superscripts 

I  =  imaginary  part 

R  =  real  part 

(  )  =  derivative  with  respect  to  time 


Introduction 

ERODYNAMICALLY-induced  vibrations  of  rotor  and 
stator  airfoils  are  one  of  the  more  common  sources  of 
high-cycle  fatigue  failure  in  gas  turbine  engines.  Destructive 
aerodynamic  forced  responses  of  fan,  compressor,  and 
turbine  blading  have  been  generated  by  a  wide  variety  of 
sources  including  upstream  blades  and/or  vanes,  distortion, 
rotating  stall,  downstream  blades  and/or  vanes,  surge,  bleeds, 
and  random  or  otherwise  unidentified  sources. 

Failure  level  vibratory  responses  occur  when  a  periodic 
aerodynamic  forcing  function,  with  frequency  equal  to  a 
natural  blade  resonant  frequency,  acts  on  a  blade  row.  The 
rotor  speeds  at  which  these  forced  responses  occur  are 
predicted  with  Campell  diagrams,  which  display  the  natural 
frequency  of  each  blade  mode  vs  rotor  speed.1  Whenever 
these  curves  cross,  aerodynamically-induced  forced  responses 
are  possible.  However,  at  the  present  time  no  accurate 
prediction  for  the  amplitude  of  the  resulting  stress  can  be 
made. 

The  prediction  of  (he  aerodynamic  forced  response 
vibratory  behavior  of  a  blade  or  vane  row  requires  a 
definition  of  the  unsteady  forcing  function  in  terms  of  its 
harmonics.  The  time-variant  aerodynamic  response  of  the 


airfoil  to  each  harmonic  of  this  forcing  function  is  then 
assumed  to  be  comprised  of  two  parts.-'  One  part  is  due  to  the 
disturbance  being  swept  past  the  non-responding  fixed  air¬ 
foils.  The  second  arises  when  the  airfoils  respond  to  this 
disturbance.  These  effects  are  modeled  by  means  of  two 
distinct  analyses.  A  linearized  small-perturbation  gust 
analysis  is  used  to  predict  the  time-variant  aerodynamics  of 
the  fixed  non-responding  airfoil  to  each  harmonic  of  the 
disturbance.  A  self-induced  unsteady  aerodynamic  analysis, 
wherein  the  airfoils  are  assumed  to  be  harmonically 
oscillating,  is  then  used  to  predict  the  additional  aerodynamic 
effect  due  to  the  airfoil  response.  These  self-induced 
aerodynamic  effects  can  be  thought  of  as  an  aerodynamic 
damping  which  can  be  either  positive  or  negative.  Reviews  of 
state-of-the-art  unsteady  aerodynamics  as  applied  to  tur¬ 
bomachines,  including  gust  and  self-induced  unsteady 
aerodynamic  analyses,  are  presented  in  Refs.  3  and  4. 

The  classical  approach  to  the  prediction  of  the  aerodynamic 
forced  response  of  an  airfoil  is  based  on  New  ton’s  second  law. 
The  gust  and  self-induced  unsteady  aerodynamic  analyses  are 
used  to  describe  the  harmonic  forces  and  moments  acting  on 
the  airfoil  in  conjunction  with  a  lumped  parameter 
description  of  the  airfoil  structural  and  inertial  properties.' 

An  alternate  approach  based  on  an  energy  balance 
technique  has  been  developed  by  Hoyniak  and  Fleeter'’  to 
predict  the  uncoupled  single-degree-of-freedom  forced 
response  of  an  airfoil.  In  this  approach,  a  balance  is 
established  between  the  energy  of  unsteady  aerodynamic  work 
and  the  energy  dissipated  by  the  airfoil. 

The  objective  of  this  paper  is  to  extend  the  energy  balance 
technique  to  include  the  more  interesting  case  of  coupled 
bending-torsion  mode  forced  response  of  an  airfoil  and,  also, 
to  demonstrate  the  effects  of  the  various  aerodynamic  and 
structural  parameters  on  this  coupled  response. 

General  Unsteady  Aerodynamic  Coefficients 

Figure  I  presents  a  schematic  representation  of  a  two- 
dimensional  airfoil  section  displaced  in  both  torsion  and 
translation  in  a  uniform  flow  with  a  superimposed  convened 
transverse  sinusoidal  gust.  Equation  ( 1 1  presents  the  complex, 
time-dependent  unsteady  lift  and  moment  per  unit  span, 
written  in  influence  coefficient  form,  for  the  gust  response 
and  the  self-induced  unsteady  aerodynamic  eases. 

L(,  =  L*  +iLl,  =  upb-'ii-'  p,,  (  q  )  j 

W,,  =  M*  +  <M\,  =  upb'i,-'  [*,.(!)] 

f-s I  =/.*  +/£.!,,  =  ir pb'^:  [/I*  -  +.4,,<<J 

Ws,  =  Mif,  +■  /A/J.,  =  irpAV’  [fl,  -  ->-  B j  ( I) 


ELASTIC  AXIS 


Fig,  I  Airfoil  and  displacement  geometry  and  notation . 


where/i  =  /i„e’-';  a  =  n„e‘-';  W=  describes  ihe  gust; 

B/,  >.  (•  ).  (.4,, .8,,)  denote  ihe  generalized  un¬ 

steady  lift  and  moment  coefficients  due  to  airfoil  translation, 
airfoil  torsion,  and  the  convected  sinusoidal  gust,  respec¬ 
tively;  u:  is  the  airfoil  natural  frequency,  and  L  is  ihe  gust 
forcing  function  frequency 

The  total  unsteady  lift  and  moment  on  the  airfoil  are  ob¬ 
tained  by  superimposing  the  gust  response  and  the  sclf- 
mduccd  aerodynamic  forces  and  moments. 

L{i)=Lt,(l)  +  Z.sl  (/) 

=  *pb’Co:  j^,;  (  —  )  +A„h„  +  A. ,nrr;je'-' 

Mil)  =Ma(l)  +WSI  (/) 

=  upb4*-’  [«,,  (  ^  )  +  B„h„  +  B.aAe"'  (2) 

It  should  be  noted  that  the  airfoil  response  occurs  at  the 
frequency  of  ihe  forcing  function.  Flcnce,  in  Eq.  (2)  the  gust 
frequency  J  has  been  utiliz.ed  in  the  specification  of  the  un¬ 
steady  aerodynamic  forces  and  moments. 


Energy  Balance 

The  equations  of  motion  describing  the  coupled  tran¬ 
slational  and  torsional  displacement  of  the  Hat  plate  airfoil 
depicted  in  Fig.  I  are  given  in  Eq.  (3). 

inh+  S„ tir  +  mu>i  U  +  ig„)h-Ht) 

S„h  +  /,,«  +  /„ ui-  ( /  +  ig„  )or  =  .W(f)  (3) 

where  m  denotes  the  mass  per  unit  span  of  the  airfoil.  S ,  is  the 
airfoil  static  moment  per  unit  span  about  the  clast ic  axis,  j,, 
and  u.,,  are  the  translational  and  torsional  mode  airfoil 
natural  frequencies  of  the  corresponding  undamped  singlc- 
dcgrce-of-frccdom  system,  respectively,  and  /.  is  the  mass 
moment  of  inertia  per  unit  span  about  the  clastic  axis.  As 
seen,  the  airfoil  structural  dynamic  system  described  in  Eq.  (3) 
is  strongly  coupled  and  the  coupling  is  associated  with  both 
the  structural  and  the  aerodynamic  characteristics  of  the 
system.  The  structural  coupling  becomes  significant  whenever 
the  aerodynamic  center  docs  not  coincide  with  the  clastic  axis, 
in  which  case  a  non-zero  value  for  the  airfoil  static  moment 
results.  The  aerodynamic  coupling  arises  because  the  self- 
induced  aerodynamic  forces  acting  on  the  airfoil  arc  a  func¬ 
tion  of  the  translational  and  torsional  motions. 

In  this  investigation,  the  coupled  system  response  is 
determined  utilizing  an  energy  balance  technique.  The  energy 
input  to  the  system  per  cycle  of  airfoil  oscillation  is  generated 
by  the  gust  and,  under  certain  conditions,  the  self-induced 
unsteady  aerodynamic  forces  amd  moments/  The  energy 
dissipation  of  the  system  per  cycle  is  associated  with  I)  system 
structural  damping,  2)  under  certain  conditions  the  self- 
induced  aerodynamic  forces  and  moments,  and  3)  the  static 
moment  term  5„.  It  should  be  noted  that,  for  the  uncoupled 
single-degree-of-freedom  case,  the  dissipation  term  5,  was 
not  considered. 

The  energy  balance  for  the  coupled  translation-torsion 
airfoil  motion  can  be  expressed  as  follows. 


(WORK)?,  -MWORK)^  +  (WORK)?  =  (WORK)?, 
(WORK)s,  +(WORK)J  -t-(WORK)"^  =  <WORK)g 


(4) 


where 


(WORK)si.  (WORK)f,  =  work  done  by  the  self-induced 
aerodynamic  forces  and  moments 
in  translation  and  torsion, 
respectively. 

(WORKlJ^ ,  (WORK)'t’u  =  work  done  by  the  airfoil  tran¬ 
slation  and  torsion  structural 
damping,  respectively. 

(WORK)*^,  iWORK);^  =work  associated  with  the  static 
moment,  S„,  in  translation  and 
torsion,  respectively. 

(WORK)J;,  (WORK)g  =  work  done  by  the  gust  aero¬ 
dynamic  forces  and  moments  in 
translation  and  torsion,  respec¬ 
tively. 

The  work  values  associated  with  the  self-induced  unsteady 
aerodynamic  forces  and  moments,  (WORK)s,  and 
(WORK)si,  have  been  written  as  dissipation  terms.  However, 
under  certain  conditions  they  may  actually  represent  energy 
input  terms.  The  sign  differentiates  between  energy  dissipa¬ 
tion  and  input.  Also,  matrix  techniques  are  utilized  to  solve 
Eq.  (4)  in  order  to  determine  the  response  of  the  airfoil 
system. 

The  unsteady  work  done  by  the  self-induced  aerodynamic 
terms,  (WORK)?,  and  ( WORK)|i ,  over  one  cycle  of  vibration 
can  be  calculated  from  Eq.  (3). 

(WORK)*,  =  <j>Re(Z.s,(r)dh| 

(WORK)si  =  <JpRelMs,  (f)dal  (5) 

where  is,  (f)  and  Ms, (r)  denote  the  unsteady  lift  and 
moment  associated  with  the  self-induced  harmonic  motion  of 
the  airfoil  (See  Eq.  (1)1 . 

Carrying  out  the  integrations  specified  by  Eqs.  I  and  5 
results  in  the  following  expressions  for  the  unsteady  work 
done  by  the  self-induced  aerodynamic  forces  and  moments. 

(WORK)?,  -  [A„)'hi,/b+  ( A„)Rhnan 

x  sin(k,  -X,,  )  -  M„  )'hyo„cos(X»  -  X„  )1 


(WORK)?,  =  tJpb4i>:  )sin(X*-X„> 

-  (fl,)1  )cos(X„  -X„)  -  ( S„  > 1  c»;, J  (6) 

The  angle  X,  denotes  the  phase  angle  between  the  incoming 
gust  and  the  resulting  translational  motion  of  the  airfoil. 
Similarly,  the  angle  X„  is  the  phase  angle  between  the  in¬ 
coming  gust  and  (he  airfoil  torsional  motion. 

The  energy  input  to  the  airfoil  is  determined  by  calculating 
(he  work  per  cycle  of  oscillation  of  the  unsteady  lift  and 
moment  associated  with  the  incoming  gust. 

(WORK)f;  =  -  |ReiLG(Mdh! 

=  -  »pbJw'hol  (Ac)'cos(X,)  -  M(,  )"sin(X, )  | 


( WORK)c  =  ^Re[,Vf<;(/)do| 


=  *pt>‘uJa0l  lBc  )'cos(X„)  -  ( Bc  )  *sin<X„  )  ( 
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Fig.  2  Comparison  of  classical  and  energy  balance  forced  response 
calculations. 
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Fig.  )  Structural  damping  effect  on  the  coupled  torsional  response 
amplitude  for  (u/w„)  near  u„. 


(7) 


The  energy  dissipated  by  the  structural  damping  over  one 
cycle  of  oscillation  can  be  calculated  as  follows. 

(WORK)^  =  ^Re(/g/,/nu>j;hd/i)  =  eg*  mu;  hi 

(  WORK)“a  =  (jp  Re  |  ig„/„u£  ador )  =  Tg„ /„ uiiuij,  (8) 


where 

An  =  [»g,mu;b)  +  |  (4*)'  r2pb,i>2  ] 

Au=  (A„  )Rr,pfr,a/lsin(X/l  ~X„ ) 

+  { (A„)[r2pbJi>2cos(\R  ~K) ) 


The  unsteady  work  associated  with  the  static  moment  S„ 
merits  some  additional  comment.  This  is  the  term  which 
couples  the  system  equations  of  motion  when  the  elastic  axis 
and  the  aerodynamic  center  do  not  coincide.  Physically,  the 
S„  term  can  be  considered  as  either  a  torsional  force  applied 
to  the  airfoil  as  the  result  of  a  unit  translation  airfoil 
displacement,  or  as  a  translational  force  resulting  from  the 
application  of  a  unit  torsional  displacement.  Thus,  for  an 
airfoil  undergoing  a  translation  mode  oscillation,  there  is  an 
apparent  torsional  load  that  can  either  do  work  on  or  extract 
work  from  the  airfoil  system.  The  phase  relationship  between 
the  translational  and  torsional  displacement  determines 
whether  this  apparent  torsional  load  adds  or  dissipates  energy 
in  the  airfoil  system. 

The  unsteady  work  per  cycle  of  oscillation  associated  with 
the  static  moment  can  be  determined  by  Eq.  (9). 

(WORK)?u  =  <|)  Re  (S„  ad/il  =S„u;a0/t<))rsin(X/, -X„) 

(WORKKo  =  ^Re[S>d„  1  =  S„i,2a0h0 sin(X„  -  X, )  (9) 

The  substitution  of  the  various  work  expressions  into  the 
energy  balance  [see  Eq.  (4)]  yields  the  following  matrix 
expression  for  the  coupled  translational  and  torsional  airfoil 
displacements. 


Aj,=  lS„i>2 t-  (fl„)Rir?ph'i)',)sin(X0-X,) 

+  |  (B/1)lr?ph4u,cos(X,  -X„) ) 

A  }3  =  Tg„  laui  +  ( Ba ) 1  rpb4u2 

Bn=  -rpb3*2  )  [  (-4c)'cos(X„)  -  (-4 G )  *sin(X* ) ) 

B2,  =  *2pb‘i,2  )'cos(X„ )- (B0  )"sm<X„ ) ) 


The  phase  angles  between  the  translational  and  torsional 
displacements,  X,  and  X„  respectively,  can  be  determined  in 
various  ways.  For  example,  they  can  be  determined  directly 
from  the  classical  solution  of  the  coupled  system  equations  of 
motion,  Eq.  (3)  (see  Ref.  5).  For  convenience  and  brevity, 
these  phase  angle  relationships  are  presented  in  Eq.  (II). 


X*  =tan 


X„  =tan 


i  (bJh,PR  —Nf,RD,  I 

[■AZ/njO*  +  Nb/D,  J 

/  fbJ„,DR-N„KD,  I 

i-N„rDr  +  N„iD,  J 


(ID 


1 

1  An 

A 12 

(*)' 

I 

B„  \ 

An 

A  22 

<0t„) 

.  1 

(10) 


where  A/*„  N„„  NhR,  N„K,  D„  and  D,  are  defined  in  terms  of 
the  airfoil  structural  properties  and  the  unsteady  aerodynamic 
forces  and  moments  in  the  Nomenclature. 


Fig.  4  Stnctaral  damping  effect  on  the  coupled  torsional  response 
for  <«/«.)  near  w*. 


I 


Fig.  S  Cascade  solidity  effect  on  the  coupled  torsional  response. 
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I'iK.  6  Cascade  solidilj  effect  on  the  cuupled  translational  response.  Hg.  (  Kffect  of  stagger  angle  on  the  coupled  translational  response. 


fig.  7  fffect  of  stagger  angle  on  the  coupled  torsional  response. 


Unsteady  Aerodynamic  Forces  and  Moments 
l„,  Af„,  L„  and  M„  denote  the  standard  form  for  the 
unsteady  aerodynamic  forces  and  moments  and  represent  the 
unsteady  lift  and  moment  in  translation  and  torsion, 
respectively,  calculated  about  the  airfoil  '/*  chord.  These 
standard  form  '/»  chord  coefficients  are  related  to  the  general 
unsteady  aerodynamic  influence  coefficients  through  the 
following  relationships. 


A*  =  Lh 

A„  =  L„-  (Zi  +  a)Lh 

Bh=Mh  -  (Zi  +a)L„ 

B„  =  M„  -  (  Zi  +a)  (/.„  +Mh)  +  (  Z,  +  a)>£* 

Bc=  Mh  (12) 

where  a  is  the  dimensionless  distance  of  the  elastic  axis 
measured  from  the  airfoil  mid-chord  (see  Fig.  1). 

Results 

To  demonstrate  the  effects  of  the  various  parameters  on  the 
coupled  forced  response  of  turbomachinery  blading,  the 
unsteady  aerodynamic  analysis  of  Ref.  7  will  be  utilized.  This 
analysis  predicts  the  gust  and  self-induced  unsteady 
aerodynamic  forces  and  moments  in  a  compressible  flow 
field.  The  structural  properties  for  this  parametric  study  are 
based  on  an  airfoil  with  a  5.08-cm  chord,  a  4%  thickness-to- 
chord  ratio,  and  an  aspect  ratio  of  3.  This  representative 
airfoil  has  a  natural  translation  mode  frequency,  <•>»,  I2°7b 
larger  than  the  natural  torsion  mode  frequency,  ua.  Hence, 
wk/w„  =  1.12.  Specific  parameters  to  be  varied  include  the 
structural  damping  value,  the  cascade  solidity,  the  interblade 
phase  angle,  the  inlet  Mach  number,  and  the  elastic  axis 
location.  In  addition,  to  clearly  demonstrate  the  coupling 
effects,  the  airfoil  stiffness  is  varied  so  as  to  make  the  tor¬ 
sional  and  translation  mode  natural  frequencies  nearly  equal 
in  value. 

To  verify  the  energy-balance  approach  for  the  prediction  of 
aerodynamically-induced  vibrations,  the  uncoupled  torsion 
mode  response  of  the  representative  airfoil  was  calculated  by 
means  of  both  the  energy  balance  technique  and  a  classical 
Newton's  second  law  approach.5  Figure  2  presents  the 
comparison  of  these  calculation  techniques.  As  shown,  the 
two  methods  yield  identical  results,  both  in  terms  of  the 
amplitude  and  the  phase  of  the  response. 
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FI*,  f  Effect  of  intcrblade  phase  angle  on  the  coupled  Ionian  made 
response. 


Fit-  10  Effect  of  Interblade  phase  angle  on  the  coupled  translation 
mode  response. 


The  coupled  forced  response  characteristics  of  the 
representative  airfoil  are  determined  utilizing  this  energy 
balance  technique.  Results  are  presented  as  a  function  of  the 
nondimensional  ratio  of  the  gust  forcing  frequency  to  the 
airfoil  natural  torsion  mode  frequency,  Hence,  values 
of  1.00  and  1.12  for  this  nondimensional  frequency  ratio 
correspond  to  the  situations  where  the  forcing  function 
frequency  is  equal  to  the  natural  torsion  and  translation  mode 
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Fig.  1 1  Inlet  Mach  number  effect  on  the  coupled  lorsionnl  response. 


frequencies,  respectively,  of  the  representative  airfoil.  Results 
presented  include  the  translational  amplitude,  the  torsion 
mode  amplitude,  and  the  phase  angles  between  the 
aerodynamic  forcing  function  and  the  resulting  airfoil 
response  as  a  function  of  the  nondimensional  frequency  ratio. 

As  anticipated,  when  the  gust  forcing  function  frequency  is 
close  to  the  natural  torsion  mode  frequency,  the  torsional 
response  is  an  order  of  magnitude  greater  than  the 
corresponding  coupled  translation  mode  response.  Thus,  only 
the  coupled  forced  torsional  response  will  be  presented  for  the 
forcing  function  close  to  the  airfoil  natural  torsional 
frequency  (w/ui„  *»  1.00),  and  only  the  coupled  forced  trans¬ 
lational  response  for  forcing  function  frequencies  near  the 
natural  translational  frequency  (&>/fa>„  ■  1.12)  will  be 
presented. 

Figures  3  and  4  represent  the  torsional  and  translational 
airfoil  responses,  respectively,  with  structural  damping  as  the 
parameter  when  the  forcing  frequency  is  close  to  u„  and  w,, 
respectively.  As  seen,  increased  structural  damping  resulting 
results  in  decreased  amplitude  of  response.  It  should  be  noted 
that  the  maximum  response  amplitudes  do  not  occur  at 
frequency  ratio  values  of  1.00  and  1.12,  even  for  the  special 
case  of  zero  structural  damping.  This  is  a  result  of  the  self- 
induced  unsteady  aerodynamic  forces  which  correspond  in 
this  case  to  aerodynamic  damping.  Thus,  even  when  there  is 
no  structural  damping  in  the  airfoil  system,  the  self-induced 
aerodynamics  generate  aerodynamic  damping. 

The  effects  of  cascade  solidity  (C/S)  and  stagger  angle  on 
the  coupled  forced  response  amplitudes  are  demonstrated  in 
Fig.  5  and  6  and  Figs.  7  and  8,  respectively.  Both  the  coupled 
torsional  and  translational  forced  response  amplitudes  in¬ 
crease  with  increasing  values  for  solidity  and  stagger  angle. 
Variations  in  solidity  result  in  a  somewhat  more  pronounced 
effect  on  the  translational  response  than  on  the  torsional 
response.  Stagger  angle  variation  has  a  greater  effect  on  the 
torsional  response  than  on  the  translational  one.  Also, 
variations  in  stagger  angle  have  a  much  greater  effect  on  the 
response  amplitudes,  both  torsional  and  translational,  than 
do  variations  in  solidity.  It  should  be  noted  that  the  design 
trend  for  modern  compressors  includes  increased  solidity  and 
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Fig.  12  Intel  Mach  number  effect  on  Ihe  coupled  translational 
response. 


stagger  angle  values.  Hence,  these  results  indicate  that  in¬ 
creased  forced  response  problems,  primarily  associated  with 
increased  stagger  angles,  might  be  anticipated. 

Figures  9  and  10  show  the  effect  of  interblade  phase  angle 
on  the  coupled  airfoil  forced  response  characteristics.  The 
effect  of  variations  in  this  parameter  on  the  coupled  torsional 
forced  response  are  quite  different  than  on  the  coupled 
translational  response.  For  the  torsion  mode,  there  is  only  an 
extreme  maximum  response  amplitude  for  a  0-deg  interblade 
phase  angle  value.  For  all  other  values  considered,  the 
response  amplitude  curves  are  relatively  Hat  and  decreased 
from  the  0  deg  case  by  approximately  75  to  85 "Ai .  The  coupled 
translational  response  amplitudes  all  exhibit  an  extreme 
maximum  response  amplitude,  i.e.  none  of  the  amplitude 
response  curves  are  relatively  flat.  Also,  although  an  in¬ 
terblade  phase  angle  value  of  0  deg  results  in  the  maximum 
forced  translational  response  amplitude,  varying  the  in- 
terblade  phase  angle  away  from  0  deg  results  in  noticabiy 
smaller  decreases  in  this  maximum  amplitude  than  were  noted 
for  the  forced  torsional  response.  In  particular,  for  both  the 
torsional  and  the  translational  cases,  the  largest  decrease  in 
the  maximum  response  amplitude  is  associated  with  changing 
the  interblade  phase  angle  value  from  0  to  -  90  deg.  However, 
the  maximum  forced  translational  response  decreased  only 
45ro.  whereas  the  corresponding  forced  torsional  response 
decreased  by  US'V  In  addition,  since  the  smallest  coupled 
torsional  and  translational  response  amplitudes  correspond  to 
negative  interblade  phase  angle  values  ( -  90  deg  in  par¬ 
ticular),  then  this  would  appear  to  be  a  desirable  forced 
response  design  condition.  In  terms  of  the  rotor-stator  in¬ 
teraction  forced  response  problem,  a  negative  interblade 
phase  angle  value  corresponds  to  a  backward  traveling  wave 
as  viewed  from  the  stator  vane  frame  of  reference  and  arises 
when  the  number  of  rotor  blades  is  greater  than  the  number  of 
stator  vanes. 

The  effect  of  the  inlet  flow  Mach  number  on  the  coupled 
torsional  and  translational  forced  response  characteristics  of 
the  representative  airfoil  is  demonstrated  in  Figs.  II  and  12, 
respectively.  Increasing  the  Mach  number,  which  also 
corresponds  to  decreasing  the  reduced  frequency  for  a  specific 
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Fig.  Id  Effect  of  elastic  axis  location  on  the  coupled  torsional 
response  for  (u/u„)  near  u„ . 
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Fig.  14  Effect  of  elastic  axis  location  on  the  coupled  translational 
response  fo r  (w/iv„ )  near  . 


gust  disturbance,  increases  both  the  coupled  torsional  and 
translational  responses  with  a  very  dramatic  increase  asso¬ 
ciated  with  an  increase  in  Mach  number  from  0.6  to  0.8.  Also, 
for  Mach  numbers  of  0.6  and  smaller,  the  increase  in  the 
maximum  response  amplitude  is  much  larger  for  the  torsional 
mode  than  for  the  translational  mode,  in  addition,  the  forcing 
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function  frequency  at  which  the  maximum  amplitude 
response  occurs  is  highly  dependent  on  the  Mach  number  for 
the  torsional  response,  but  is  almost  independent  of  this 
forcing  frequency  for  the  translational  response  case. 

The  effect  of  moving  the  elastic  axis  location  to  the  '/<- 
chord  and  the  ’/.-chord  locations  on  the  coupled  forced 
response  of  the  representative  airfoil  is  presented  in  Figs.  13 
through  16.  These  results  indicate  that  not  only  the  response 
amplitude  but  also  the  forcing  function  frequency  at  which 
ihe  maximum  amplitude  response  occurs  is  strongly 
dependent  on  the  position  of  the  elastic  axis.  The  maximum 
torsion  and  translation  response  amplitudes  are  found  with 
the  elastic  axis  at  the  ’/.-chord  location.  Shifting  the  elastic 
axis  forward  to  the  '/.-chord  location  results  in  significantly 
decreasing  both  the  torsional  and  translational  response 
amplitudes  near  the  airfoil  torsional  natural  frequency,  but 
only  slightly  decreasing  these  response  amplitudes  near  the 
airfoil  translational  natural  frequency.  Also,  shifting  the 
elastic  axis  location  aft  from  mid-chord  to  ’/.-chord  has  the 
result  of  making  the  translational  response  in  the  neigh¬ 
borhood  of  the  natural  torsional  frequency  (w/w„  »  l  .0)  the 
same  order  of  magnitude  as  the  torsional  response  amplitude, 
indicating  the  significant  additional  coupling  that  arises  when 
the  elastic  axis  and  the  airfoil  center  of  gravity  do  not  coin¬ 
cide.  Similar  results  are  noted  when  the  forcing  function 
frequency  is  close  to  the  natural  translational  frequency. 

To  demonstrate  the  interaction  of  two  vety  closely-spaced 
torsional  and  translation  modes,  the  bending  stiffness  of  the 
representative  airfoil  was  altered  to  make  the  ratio  of  tran¬ 
slation  to  torsion  natural  frequencies  1.02.  The  variation  with 
interblade  phase  angle  values  of  these  closely  coupled  forced 
torsional  and  translational  responses  was  then  investigated. 
The  results  are  presented  in  Figs.  1 7  and  18.  The  coupled 
torsional  response  results  (Fig.  17}  show  that,  at  a  frequency 
ratio  of  0.990,  the  response  decreases  as  the  interblade  phase 
angle  is  varied  from  0  to  +90,  -90,  and  180  deg.  At  a 
frequency  ratio  of  1 .01 5  the  response  increases  by  up  to  f>Vrt . 
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Fig.  17  Inttrblide  phase  ingle  effect  on  the  coupled  torsional 
response  for  nearly  equal  natural  frequencies. 


Figure  18  presents  the  variation  of  the  forced  translational 
response  with  the  interblade  phase  angle  as  the  parameter.  In 
contrast  to  the  torsional  response  results,  minimal  trans¬ 
lational  response  is  noted  when  the  forcing  function 
frequency  is  near  the  airfoil  natural  torsional  frequency. 
However,  when  the  forcing  function  frequency  is  near  the 
airfoil  natural  translational  frequency  (£>/«„  ■  1.015),  the 
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forced  translational  amplitude  is  significant  and  increases  as 
the  interblade  phase  angle  value  is  varied  from  0  to  180  deg.  A 
comparison  of  Figs.  17  and  18  at  a  frequency  ratio  close  to 
1.015  reveals  a  correspondence  between  the  forced  torsional 
response  and  the  translational  one.  This  is  caused  by  the  high 
degree  of  coupling  that  exists  between  the  translational  and 
torsional  modes  for  the  altered  airfoil. 

Summary  and  Conclusions 

In  this  investigation,  the  coupled  structural  dynamic  airfoil 
system  forced  response  was  determined  utilizing  an  energy 
balance  technique.  The  energy  input  to  the  system  per  cycle  of 
oscillation  was  generated  by  the  gust  and,  under  certain 
conditions,  the  self-induced  aerodynamic  forces  and 
moments.  The  energy  dissipation  per  cycle  of  oscillation  was 
associated  with  system  structural  damping,  under  some 
conditions  the  self-induced  aerodynamic  forces  and  moments, 
and  the  static  moment  term  5„.  It  should  be  noted  that,  for 
the  uncoupled  single-degree-of-freedom  case,  the  dissipative 
coupling  mechanism,  S„,  is  not  considered. 
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The  effects  of  the  various  aerodynamic  parameters  on  the 
coupled  translational  and  torsional  mode  forced  response  of  a 
representative  airfoil  were  then  considered.  The  study  showed 
the  increased  coupling  between  the  torsion  and  translation 
modes  as  the  corresponding  undamped  natural  frequencies 
approach  one  another.  It  was  also  demonstrated  that  the 
coupled  torsional  and  translational  forced  response  am¬ 
plitudes  of  a  representative  airfoil  increased  with:  I )  decreased 
structural  damping,  2)  increased  solidity  values,  3)  increased 
stagger  angles,  4)  increased  inlet  flow  Mach  numbers,  5) 
interblade  phase  angle  values  corresponding  to  forward 
traveling  waves  for  the  rotor-stator  interaction  case,  and  6) 
shifting  of  the  elastic  axis  location  aft. 

It  should  be  noted  that  variations  in  the  above  parameters 
did  not  affect  the  magnitude  of  the  resulting  forced  response 
equally  nor  did  they  always  have  equivalent  effects  on  the 
torsional  and  translational  response  modes.  For  example, 
increasing  the  stagger  angle  and  shifting  the  elastic  axis  aft 
resulted  in  significantly  larger  increases  in  all  of  the  response 
amplitudes  than  did  increasing  the  solidity.  Also,  the  in- 
terbiade  phase  angle,  for  example,  had  a  much  greater  effect 
on  the  coupled  torsion  mode  responses  than  on  the  trans¬ 
lational  ones.  In  addition,  the  forcing  function  frequencies  at 
which  the  maximum  torsional  and  translational  responses 
occurred  did  not  generally  correspond  to  the  airfoil  natural 
torsional  or  tanslational  frequencies.  This  was  because  the 
self-induced  aerodynamic  forces  and  moments  result  in  an 
aerodynamic  damping  effect  analogous  to  structural  damp¬ 
ing. 

Acknowledgment 

The  Air  Force  Office  of  Scientific  Research  and  Capt. 
Michael  Francis,  who  served  as  the  project’s  technical 
monitor,  are  gratefuly  acknowledged  for  their  support. 


References 

'••The  Aerothermodynamics  of  Aircraft  Gas  Turbine  Engines,” 
AFAPL-TR  78-52.  July  1978. 

JFleeter,  S..  Jay,  R.L.,  and  Bennett,  W.A.,  “Rotor  Wake 
Generated  Unsteady  Aerodynamic  Response  of  a  Compressor 
Stator,”  ASME  Journal  of  Engineering  for  Power,  Vol.  100.  No.  4, 
Oct.  1978,  pp.  664-675. 

’Fleeter,  S.,  "Aeroelasticity  Research  for  Turbomachine  Ap¬ 
plication,"  Journal  of  Aircraft,  Vol.  16,  May  1979,  pp.  320-326. 

J  Platzer,  M.F.,  “Unsteady  Flows  in  Turbomachines— A  Review  of 
Recent  Developments,"  AGARD-CP-227,  Sept.  1977. 

’Fung,  Y.C.,  An  Introduction  to  the  Theory  of  Aeroelasticity, 
Dover  Publications,  New  York,  N.Y.,  1969. 

6Floyniak,  D.  and  Fleeter,  S.,  "Prediction  of  Aerodynamically 
Induced  Vibrations  in  Turbomachinery  Blading,"  ASME  Symposium 
in  Fluid/Slruclure  Interaction  in  Turbomachinery,  Nov.  1981 . 

’Fleeter,  S.,  “The  Fluctuating  Lift  and  Moment  Coefficients  for 
Cascaded  Airfoils  in  a  Nonuniform  Compressible  Flow,”  Journal  of 
Aircraft.  Vol.  10,  Feb.  1973,  pp.  93-98. 


Ji-86 


DTIC 


